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Abstract 

Let K = (Ki, . . . , K n ) be an n-tuple of convex compact subsets in the Euclidean space 
R™, and let V(-) be the Euclidean volume in R™. The Minkowski polynomial Vk is defined 
as Vk(Ai, . . . , A„) = V(\\Ki+ 7 • • • , +X n K n ) and the mixed volume V{K\, . . . , K n ) as 

V(K u ...,K n ) = ^ Vk(Ai^i+,---,+A K X»). 
aXi...oX n 

Our main result is a poly-time algorithm which approximates V(Ki, . . . ,K n ) with multi- 
plicative error e™ and with better rates if the affine dimensions of most of the sets Ki are 
small. Our approach is based on a particular approximation of log(V(Ki, . . . ,K n )) by a 
solution of some convex minimization problem. We prove the mixed volume analogues of 
the Van der Waerden and Schrijver- Valiant conjectures on the permanent. These results, 
interesting on their own, allow us to justify the abovementioned approximation by a convex 
minimization, which is solved using the ellipsoid method and a randomized poly-time time 
algorithm for the approximation of the volume of a convex set. 
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1 Introduction 



Let K = (Ki, . . . , K n ) be a n-tuple of convex compact subsets in the Euclidean space R n , 
and let V(-) be the Euclidean volume in R n . It is a well known result of Hermann Minkowski 
(see for instance [5]), that the value of Vk (Aii^i + • • • X n K n ), where "+" denotes Minkowski 
sum, and XK denotes the dilatation of K with coefficient A, is a homogeneous polynomial of 
degree n in nonnegative variables Ai...A n (called the Minkowski polynomial). The coefficient 
V(K±, . . . , K n ) of Ai • A2 • • • • A n is called the mixed volume of K\...K n . Alternatively, 

V{K U ...,K n ) = —JZL—VKfaK! + • • • + \ n K n ). 

OA\...0\ n 

Mixed volume is known to be monotone [5], namely Ki C Lj, for i = 1, n, implies V(K%, . . . , K n ) < 
V{L\...L n ). In particular, it is always nonnegative and therefore all the coefficients of the 
Minkowski polynomial Vk are nonnegative real numbers. 

The corresponding Brunn-Minkowski theory, which is the backbone of convex geometry 
and its numerous applications, is about various implications of the fact that the functional 
(Vk(Ai.Ki + ■ ■ ■ \ n K n ))n is concave on the nonnegative orthant i?" = {(Ai,...A n ) : Aj > 0}. 
Its generalization, Alexandrov-Fenchel theory, is based on the fact that the functionals 
( 9Xl dXk Vk(0, ..,0, Afc+i, A n )"- fe are concave on R^T for all 1 < k < n — 1. 

The problem of computing the mixed volume of convex bodies is also important for combi- 
natorics and algebraic geometry [7]. For instance, the number of toric solutions to a system of 
n polynomial equations on C n is upper bounded by — and for a generic system, equal to — the 
mixed volume of the Newton polytopes of the the corresponding polynomials. This remarkable 
result, called the BKK Theorem, is covered, for instance, in [25] and [5]). 



1.1 Previous Work 

The BKK Theorem created an "industry" of computing (exactly) the mixed volume of integer 
polytopes and its various generalizations; most algorithms in the area are of exponential runing 
time ([H],[9], [H] and many more). Most researchers in the "industry" don't bother to formally 
write down the complexity, rather they describe the actual amount of the computer time. 
Although there was a substantial algorithmic activity on the mixed volume of polytopes prior 
to [7], the paper [7J was the first, to our knowledge, systematic complexity-theoretic study in the 
area. It followed naturally upon the famous FPRAS algorithms [6] for volumes of convex bodies, 
solved several natural problems and posed many important hard questions. The existence of 
FPRAS for the mixed volume even for polytopes or ellipsoids is still an open problem. 

Efficient polynomial-time probabilistic algorithms that approximate the mixed volume ex- 
tremely tightly (within a (1+e) factor) were developed for some classes of well-presented con- 
vex bodies [7J. The algorithms in [7J are based on the multivariate polynomial interpolation 
and work if and only if the number k of distinct convex sets in the tuple K is "small", i.e. 
k = 0(log(n)). 

The first efficient probabilistic algorithm that provides a n°( n )-factor approximation for arbi- 
trary well-presented proper convex bodies was obtained by Barvinok [2]. Barvinok's algorithms 
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start by replacing convex bodies with ellipsoids. This first step already gives n°( n )-factor in the 
worst case. After that the mixed volume of ellipsoids is approximated with a simply exponential 
factor c n by two randomized algorithms; one of those deals with approximation of the mixed 
discriminant. 

The question of existence of an efficient deterministic algorithm for approximating the mixed 
volume of arbitrary well-presented proper convex bodies with an error depending only on the 
dimension was posed in [7J. The authors quote a lower bound (Barany and Fiiredi bound) 

n 

PQ of (p, 2 for the approximation factor of such an algorithm. (Notice that Barvinok's 

randomized algorithm [2 J does not beat the Barany and Fiiredi bound.) 

Deterministic polynomial-time algorithms that approximate the mixed volume with a factor 
of were given, for a fixed number of distinct proper convex bodies in K = {K%, . . . , K n ), 

i n [2] j [Z] ■ Finally, a deterministic polynomial-time algorithm that approximates the mixed vol- 
ume with a factor of 

n O(n) in 

the general case of well-presented compact convex sets was given 
in [13], [2]. Let A = (A%, ...,A n ) be an n-tuple ofnxn complex matrices; the corresponding 
determinantal polynomial is defined as Detx{\\, A n ) = det(^i<i< ri , AjAj). The mixed dis- 
criminant is defined as D(A\, A n ) = dx ® dx Bet a(Ai, X n ). Similarly to the randomized 
algorithm from [2] , the algorithm in [T3J , [13] reduced the approximation of the mixed volume of 
well-presented compact convex sets to the approximation of the mixed volume of ellipsoids; this 
first step gives an factor in the worst case. And the mixed volume of ellipsoids is approx- 

imated by (D(A\, A n )) 2 of the corresponding positive semidefinite matrices A4 >z 0. This 
second step adds y/3 to the multiplicative approximation error (see inequality ()16|) below). 

The approximation of the mixed discriminant has also been relaxed to a convex optimization 
problem (geometric programming). In order to prove the accuracy of the convex relaxation, 
the author proved in [15] the mixed discriminant analogue of the Van der Waerden conjecture 
on permanents of doubly stochastic matrices [21], which was posed by R. V. Bapat in [3j. 

To summarize, the interpolational approach from [7] is limited by the restriction that the 
number of distinct convex sets is 0(log(n)); previous approaches [2], [13], p3] can't give the 
simply exponential approximation factor c n because of the initial approximation of convex sets 
by ellipsoids. 



1.2 Our Approach 

Assume, modulo deterministic poly-time preprocessing [7J, that the mixed volume V(K\, . . . , K n ) > 
0. We define the capacity of the volume polynomial Vk as 

Ca P (V K )= inf •■•'*") (1) 

x t >0:l<i<n !!!<;<„ 

Since the coefficients of the volume polynomial Vk are nonnegative real numbers we get the 
upper bound 

V(Ki,...,K n ) < Ca P (V K ). 
The trick is that log(Cap(Vk)) is a solution of the following convex minimization problem 

log(Cap(V K )) = inf log(Vk(e^,...,e^)). (2) 
yi + ...+y n =0 
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Recall that the functional log(p(e yi , e v ")) is convex on R n if p(x\, x n ) is any polynomial 
with nonnegative coefficients. More generally, a sum of log- convex functionals is log-convex. 

We view Cap(Vk) as an approximation for the mixed volume V{K\, . . . ,K n ). To justify 
this we prove the lower bound 

71 ' 

V{K X , ...,K n )> —Cap(V K ) « e- n Cap(V K ), 
n n 

which is the mixed volume analogue of the Van der Waerden conjecture on the permanent. 
We also present better upper bounds when "most" of the convex sets K{ have small afflne 
dimension, which are analogues of the Schrijver- Valiant conjecture, posed in [26] and proved in 
[27]. 

The idea of our approach is very similar to our treatment of H-Stable polynomials in 
|17j . Recall that a homogeneous polynomial p(x%, ...,x n ) with nonnegative coefficients is called 
H-Stable if p(z±, z Tl ) ^ provided that the real parts Re(zi) > 0, 1 < i < n. 

Not all Minkowski polynomials Vk are H-Stable: any univariate polynomial with nonnega- 
tive coefficients S(x) = Y^,o<i<n CD a « xl such that af > aj_iai_i, 1 < i < n — 1 can be presented 
as S(x) = V(A+xB) for some convex compact subsets (simplexes) A, B C R n [28J. Fortunately, 
a modification of the inductive proof in [T7j works for Minkowski polynomials and presented in 
the next Section. 

After establishing the mixed volume analogues of the Van der Waerden and the Schrijver- 
Valiant (permanental) conjectures, we present a randomized poly-time algorithm to solve the 
problem ([2]) based on the ellipsoid method and randomized poly-time algorithms for the vol- 
ume approximation. Together with the mixed volume analogues of the Van der Waerden 
conjecture this gives a randomized poly-time algorithm to approximate the mixed volume 
V(K\, . . . ,K n ) within relative accuracy e n . Notice that, in view of the Barany and Fiiredi 
bound, this cannot be achieved by a deterministic poly-time oracle algorithm. We use the 
ellipsoid method because of its robustness: we deal essentially with a random oracle which 
computes log (Vk(e yi , e Vn )) with an additive small error e ; we use this oracle to get an 
approximation of the gradient of log (Vk^ 1 , ...,e Vn )). 

2 Van der Waerden and Schrijver- Valiant conjectures for the 
mixed volume 

Consult Appendix A for the proofs of results in this section. 

2.1 The mixed volume analogue of the Van der Waerden- Falikman-Egorychev 
inequality 

Theorem 2.1: 
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1. Let K = (Ki, . . . , K n ) be a n-tuple of convex compact subsets in the Euclidean space R n . 
The mixed volume V(K\, . . . , K n ) satisfies the next lower bound: 

n ' 

V(K 1 ,...,K n )>—Cap(V K ) (3) 

n n 

2. The equality in ([3|) is attained if and only if either the mixed volume V{K\, . . . , K n ) = 
or Ki = diKi + bi : ai > 0,bi e R n ; 2 < i < n. 

2.2 The mixed volume analogue of the Schrijver- Valiant conjecture 
Definition 2.2: 

1. Let n > k > 1 be two integers. We define the univariate polynomial 



sv n<k (x) = 1 + 

Ki<k 
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Note that sv n>n (x) = (l + ^) n . Next we define the functions : 

X(n,k) = (unn( S -^l)Y 1 (4) 

\x>0 x ) 



Remark 2.3: It was observed in |16| that 



X(k,k)=g(k)=: (^) fc \fc>l; II 9(k) = ^ (5) 

l<k<n ' 

The following inequalities are easily verified: 

A(n, k) < A(n, I) : n > k > I > 1; X(m, k) > X(n, k) : n > m > k. (6) 
It follows that 

/ Y, <•< t-\ 1 

X(oo, k) =: lim X(n,k) = I mm =^= — —I (7) 

-i 



The equality A(n, 2) = (1 + v^V^) > (1 + v^)^ 1 follows from basic calculus. I 



2. Let n > m > 1 be two integers. An univariate polynomial with nonnegative coefficients 
R(t) = J2o<i<m a it l 1S called n — Newton if it satisfies the following inequalities : 

(The standard Newton's inequalities correspond to the case n = m are satisfied if all the 
roots of p are real.) 
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The main mathematical result in this paper is the following theorem: 

Theorem 2.4: Let K = (K\, ■ ■ ■ ,K n ) be a n-tuple of convex compact subsets in the Euclidean 
space R n and aff(i) be the affine dimension of Ki, 1 < i < n. Then the following inequality 
holds: 

Cap(V K )>V(K 1 ,...,K n )> J] X(i,D(i))Cap(V K );D(i)=mm(i,aff(i)),l<i<n. (9) 

l<i<n 

Corollary 2.5: Suppose that aff(i) < k : k + 1 < i < n then 

V(K 1: ...,K n )> pA(n, k) n - k Cap(V K ). (10) 
If k = 2 we get the inequality 

V(Kx, . . . ,K n ) > ~(1 + V2) 2 - n Cap(V K ) (11) 

2.3 Comparison with previous results 

The inequality @ is an analogue of the famous Van der Waerden conjecture [21], proved in 
[10] and [8], on the permanent of doubly-stochastic matrices. Indeed, consider the "boxes" 
Ki = {(xi, x n ) : < Xj < A(i,j),l < j < n}. Then the mixed volume is equal to the 
permanent, 

V(K) =V(K 1 , . . . ,K n ) = Perm(A), 

and if the n x n matrix A is doubly-stochastic then CapiVn) = 1. 

Though this mixed volume representation of the permanent has been known since the pub- 
lication of [S] if not earlier, the author is not aware of any attempts prior to the present paper 
to generalize the Van der Waerden conjecture to the mixed volume. We think that our version, 
stated in terms of the capacity, is most natural and useful. 

The inequality (fTO]) is an analogue of Schrijver's lower bound [27], [TTJ on the number of 
perfect matchings in /c-regular bipartite graphs: affine dimensions play role of the degrees of 
vertices. 

The reader familiar with [TTJ can recognize the similarity between inequalities ([9]), ([3]), (|10p 
and the corresponding inequalities in [TTJ, proved for H-Stable polynomials. The method of 
proof in the present paper is also similar to the one in [TTJ (in spite of the fact that not all 
Minkowski polynomials Vk are H-Stable). But we get worse constants: for instance, if k = 2, 
in the notation of (|10p . then in the H-Stable case we get the factor 2~ n+1 instead of the 
|(1 + V2) 2 ~ n obtained in this paper. Whether the latter factor is asymptotically sharp is an 
open problem. 
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2.4 The idea of our proof 



Our proof of Bapat's conjecture in [15], i.e. of Van der Waerden conjecture for the mixed 
discriminant, is an adaptation of Egorychev's proof in [8]. In contrast, the proofs in [T7] and 
in the present paper have practically nothing in common either with Egorychev's proof or with 
Falikman's proof in |10j . 



1. How do we prove the lower bounds? 

We follow the general approach by the present author, introduced in |17j . 
We associate with the Minkowski polynomial Vk a sequence of polynomials: 

Qn—i 

Qn = V K ,qi(xi, ...,Xi) = — q n (xi, ...,Xi,0, ...,0), 1 < % < n-1. 

dx i+1 ...dx n 

Note that qi(x) = V(K\, . . . ,K n )x. Everything follows from the next inequality 

Cap(qi) > A (i + l,min(i + 1, aff(i + 1))) Cap(q i+ i), 1 < i < n— 1. 

Not surprisingly, we do use the (still hard to prove) Alexandrov-Fenchel inequalities to 
prove this crucial inequality. (In contrast, the H- Stable case in [17] required just the 
elementary AG inequality). 

2. How do we prove the uniqueness? 

The uniqueness proofs in [8j and |15j are critically based on the known characterization of 
the equality cases in the Alexandrov inequalities for the mixed discriminant. In the case 
of the Alexandrov-Fenchel inequalities for the mixed volume such a characterization is not 
known. Luckily, the method of our proof of the lower bound ([3]) allows us to use the well 
known characterization of equality in the (much simpler) Brunn-Minkowski inequality. 
The uniqueness proof in the present paper is very similar to the uniqueness proof in |17j . 
The fundamental tool in [T7] was Garding's famous (and not hard to prove) result on the 
convexity of the hyperbolic cone. 



3 Convex Optimization Relaxation of the Mixed Volume 

Inequalities ([U El [TUl) justify the following strategy for approximation of the mixed volume 
V(K\, . . . , K n ) within a simply exponential multiplicative factor: solve the convex optimization 
problem ([2]) with an additive O(l) error. We follow here the approach from [13|. 1 14] which dealt 
with the following problem: 

\og(Cap(Det A )) = inf log [ det( V e Vi Ai) I : A { y 0. (12) 
vi+...+j/n=o \ frt ) 

\ l<j<n ) 

The main difference between the two problems is that the value and the gradient of determi- 
nantal polynomials can be exactly evaluated in deterministic polynomial time. The case of the 
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Minkowski polynomials Vk requires some extra care. Yet, this is done using standard and well 
known tools. This section of the paper is fairly routine and can be easily reproduced by any 
convex optimization professional. 

We now give an overview of the main points: 



3.1 A brief overview 

1. Representations of convex sets and a priori ball for the convex relaxation: we 

deal in this paper with two types of representations. First, similar to [7] we consider 
well-presented convex compact sets; second, motivated by algebraic applications and the 
BKK theorem, we consider integer polytopes given as a list of extreme points. In both 
cases, we start with deterministic poly-time preprocessing which transforms the initial 
tuple K G R n into a collection of indecomposable tuples 

KW eR dl ,...,K® eR^; ^ d j = n 

l<j<i 

such that the mixed volume V(K) = IIi<j<i V(K^). The tuple K is indecomposable if 
and only if the minimum in ([2]) is attained and unique. The preprocessing is essentially 
the same as in [14) . 

After this preprocessing we deal only with the indecomposable case and get a priori ball 
which is guaranteed to contain the unique minimizer of ([2]). The radius of this ball is 
expressed in terms of the complexity of the corresponding representation: 
r < 0(n 2 (log(n)+ < K >)), where < K > is the complexity of the initial tuple K. 

This part is fairly similar to the analogous problem for (|12p treated in |14j . 

2. Lipschitz Property and Rounding: In the course of our algorithm we need to evaluate 
the volumes V(e Vl Ki + ... + e Vn K n ) and the mixed volumes V(Ki, B, B), B = e Vl K\ + 
. . . + e Vn K n . This requires a well-presentation of the Minkowski sum B = e Vl K\ + ... + 
e y "K n . Given the well-presentation of K one gets a well-presentation of aiKi + ...+a n K n if 
the sizes of positive rational numbers en are bounded by poly(n, < K >) [7j. Therefore we 
need a rounding procedure, which requires us to keep only a "small" number of fractional 
bits of yi (integer bits are taken care of by a priori ball). This can be done using the 
Lipschitz property (fT8l) of log (V(e Vl Ki + . . . + e Vn K n )) (which is just Euler's identity for 
homogeneous functionals) and its partial derivatives, proved in Lemma (l3.3p . 

3. Complexity of our Algorithm: We give an upper bound on the number of calls to 
oracles for Minkowski sums a\K\-\- ...+a n K n . The number of calls to oracles for the initial 
tuple K will be larger but still polynomial (see the discussion in [4] in the context of surface 
area computation, which is, up to a constant, the mixed volume V(Ball n (l), A, A)). 

4. Ellipsoid Method with noisy first order oracle: let g(.) be a differentiable convex 
functional defined on the closed ball Ball n (r) = {X € R n :< X,X >< r 2 } and Var(g) = 
m&X-xeBall n (r) ffC^O ~~ m ^ n XGBaii n (r) ffO^O- The standard version of the ellipsoid method 
requires exact values of the function and its gradient. Fortunately, there exists a noisy 
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version [22], which needs approximations of the value g(X) and of the gradient (\jg)(X) 
such that 

sup YtXeBaUn(r) mX)+ < &T)(X),Y-X >)-(g(X) + ( V g)(X),Y-X >)| < 5Var(g). 

In our case, g(yi, y n ) = log(Vk(e yi , e Vn )). We get the additive approximation of 
g(yi, y n ) using a FPRAS for the volume approximation and the additive approxima- 
tion of (xjg)(X) using FPRAS from [7] for approximating the "simple" mixed volume 
(generalized surface area) V{Ki, K, K),K = J2i<i< n e Vi Ki- 

3.2 Representations of convex compact sets 

Following [7] we consider the following well-presentation of convex compact set C R n , 1 < 
i < n: a weak membership oracle for K together with a rational nx n matrix A{ and a rational 
vector Yi G R n such that 

Yi + Ai{Ball n {\)) cKiCYi + nV : n~+lA i (Ball n (l)) (13) 

We define the size < K > as the maximum of bit sizes of entries of matrices Ai, 1 < i < n. 
Since the mixed volume V{K\, . . . , K n ) = V{K\+{—Y\}, K n +{— Y n }, we will assume WLOG 
that Yi = 0, 1 < i < n. This assumption implies the following identity for affine dimensions 

^) = Rank(J2 AAj), S C {1, .., n}. (14) 

ies ies 



Definition 3.1: An n-tuple K = {K\ , . . . , K n ) of convex compact subsets in R n is called 
indecomposable if aff(J2ieS ^i) > Card(S) : S C {1, n}, 1 < Card(S) < n. 
We consider, similar to [2], n(n — 1) auxiliary n-tuples K y , where K y is obtained from K by 
substituting K{ instead of Kj . Notice that 

V(x 1 K 1 + ... + x n K n ) = x 1 X2...x n (V(^) + \ E ^V(K ij )) + - (15) 

l<i,j<n 3 

I 

It follows from (114p that the n-tuple K = {K\, . . . ,K n ) of well-presented convex sets is inde- 
composable iff the n-tuple of positive semidefinite matrices Q = (Qi...Q n ) '■ Qi = MAj is fully 
indecomposable as defined in |14j . which implies that indecomposability of K is equivalent to 
the inequalities V(K lJ ) > : 1 < i, j < n. Here V^K 13 ) stands for the mixed volume of the 
n-tuple K 1 - 5 . 



It was proved in \1J$ that an n-tuple of positive semidefinite matrices Q = (Qi-.-Q n ) is 
indecomposable if and only if there exists an unique minimum in the optimization problem 

inf log(det( E e*Qi)). 



s 



In the same way, an n-tuple K = (K\, . . . ,K n ) of convex compact subsets in R n is indecom- 
posable if and only if there exists an unique minimum in the optimization problem |1|). 

Applying the decomposition algorithm from Section 2 in [H] to n-tuple of positive semidef- 
inite matrices Q = (Qi--.Q n ), we can, by deterministic poly-time preproprocessing, determine 
whether or not the n-tuple of convex compact subsets K is indecomposable and if not, factor the 
mixed volume as V(K) = Yli<j< m < n ^C^-j)- Here the n(j)-tuple Kj = (Kj t i, Kj >n ^) C R n ^ 
is well presented and indecomposable, J2i<j<m n (i) = n an d the sizes < Kj > < < K > 
+poly(n). 



Based on the above remarks, we will deal from now on only with indecomposable well- 
presented tuples of convex compact sets. Moreover, to simplify the exposition, we assume 
WLOG that the matrices A{ in (|13p are integer. 



Let &a be the ellipsoid A{Ball n {\)) in R n . The following inequality, proved in [2], connects 
the mixed volume of ellipsoids and the corresponding mixed discriminant: 

3- s Pv n D*{A 1 {A 1 ) T ,...,A n (A n ) T ) <V(£ Al ...£ An ) < v n D*{A 1 {A 1 ) T ,...,A n (A n f). (16) 
Here v n is the volume of the unit ball in R n . 



3.3 Properties of volume polynomials: Lipschitz, bound on the second deriva- 
tive, a priori ball 

Proposition 3.2: 

1. Lipschitz Property. 

Let p(xi, be a nonzero homogeneous polynomial of degree n with nonnegative coef- 

ficients, Xi = e Vi . Then 

® i i t v-i v n \\ •••' x «) e?/! 

— log(p(e yi ,...,e Vn )) = —^- r (17) 

% p(x 1 ,...,x n ) 

It follows from the Euler's identity that J2i<i< n ~§^~ \og{p(e Vl , ...,e Vn )) = n, therefore the 
functional f(yi,...,y n ) = log(p(e yi , e Vn )) is Lipschitz on R n : 

+ 5 1 ,...,y n +S n ) - f(yi,...,y n )\ < ™||A||oo < n||A|| 2 (18) 

2. Upper bound on second derivatives. 

Let us fix real numbers yi, —,yi-i,yi+i, —,y n an d define univariate function q(yi) = 
log(y K (e yi ,-,e 2/ S-,e 2/ "))- Notice that e*M = Eo<j<aff( Ki ) a, > 0. 

Proposition 3.3: 

0<q"(y)<aff(K l ) . (19) 
(Lemma W.3\) in Appendix B proves a more general inequality.) 
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3. A Priori Ball result from [14] . 

Let p G Hom + (n,n), p(xi, ...,x n ) = Xix 2 ...x n (a + \ T,i<i^j< n &lJ fj) + •••• 

Assume that minKijtj^b 1 ^ =: Stf(p) > 0. Then there exists an unique minimizer 

{zi, ...,z n ),J2i<i<n z i =° such that 

\og(p(zi,...,z n )) = min (log(p(e ?/1 ,...,e ?/n )) . 

Ei< ! <„^=° 

Moreover, 

The next proposition directly adapts Lemma 4.1 from [14J to the mixed volume situation, using 
Barvinok's inequality (|16f) . 

Proposition 3.4: Consider an indecomposable n-tuple of convex compact sets K = (Ki, . . . , K n ) 
with the well-presentation Ai(Ball n (l)) C Ki C y + ri\fn + lAi(Ball n (l)), 1 < i < n with inte- 
ger n x n matrices A^. Then the minimum in the convex optimization problem ([JjJ is attained 
and unique. The unique minimizing vector (zi, z n ), J2i<i< n z i = *-* satisfies the following 
inequalities: 

\zi - zj\ < 0(nt(log(n)+ < K >)); \\ Zl , ...,z n )\\ 2 < 0(n 2 (log(n)+ < K >)) (21) 
In other words the convex optimization problem (0|) can be solved on the following ball in R n ~ l : 
Apr(K)={(z 1 ,...,z n ):\\z 1 ,...,z n )\\ 2 <0(n 2 (log(n)+<K>)), £ Zi = 0} (22) 

l<i<n 

The following inequality follows from the Lipschitz property U8\): 

\log(V K (e y \...,e yn ) -log(V K (e ll ,...,e ln )\ < 0(n 3 (log(n)+ < K >)) : Y,L e Apr(K). (23) 

3.4 Ellipsoid method with noisy first order oracles 

We recall the following fundamental result [22] : 



Let f(Y) be differ entiable convex functional defined on the ball Ball n {r) = {Y £ R n :< 
Y, Y >< r 2 } of radius r. Let Var(f) = max yeMln(r) f(Y) - miny gjBa ; /n(r ) f(Y). Assume that 
at each vectorY £ Ball n {r) we have an oracle evaluating a value g(Y) such that \g(Y) — f(Y)\ < 
0.25Var(f) and the vector gr(Y) G R n such that \\gr(Y) - (v/)(X)l|2 < 0.26r~ 1 Var(f) (here 
(V f)(Y) is the gradient of f evaluated at Y). Then the Ellipsoid method finds a vector 
Z £ Ball n {r) such that f{Z) < mm Y& Baii n (r) f(Y) + eVar(f),e > 5. The method requires 
0(n 2 log(^^)) oracle calls plus 0{n 2 ) elementary operations to run the algorithm itself. 
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3.5 Putting things together 



Here we take advantage of randomized algorithms which can evaluate log(yoZ(i^)),for a well- 
presented convex set K, with an additive error e and failure probability 5 in 0(e~ k n l log(|)) 
oracle calls. For instance, the best current algorithm [20] gives k = 2, 1 = 4. We will need below 
to evaluate volumes V(J2i<i< n x %^i)- in our case the functional / = log(VK(e yi , e Vn ) defined 
on ball Apr(K) of radius 0(n 2 (log < K >)) with the variance Var(f) < 0(n 3 (log(n)+ < 
K >)). Theorem 12.11 gives the bound: 

log(V(K))<( min /(F)) < log(F(K)) + log(— ) « log(V(K)) + n. (24) 
YeApr(K) n\ 

Therefore, to approximate the mixed volume V(K) up to a multiplicative factor of e n it is 
sufficient to find Z £ Apr(K) such that f(Z) < mm Y( zA P r(K) /00 + 0(1). In order to get 
that via the Ellipsoid method we need to approximate \og(Viz(e Vl , ...,e Vn )) with the additive 
error OiVar(f)^ 1 ) = 0(n~ 3 (log(n)+ < K >) _1 ) and its gradient with the additive 1% error 
0(n- 2 (log(ra)+ < K >)~ x ). 



1. Approximation of log(Vk(e i/1 , e Vn )) with failure probability 5. The complexity 
is O(n 10 (log(n)+ < K >) 2 log(5" 1 )) 

2. Approximation of the partial derivatives. Let Xj = e Vi and recall that the partial 
derivatives are 

oyi Vk{xi, ...,x n ) 

Suppose that < 1 — a < ^ < 1 + a. It follows from the Euler's identity that Y^i<i< n 

Pi\ < a. If a = 0(n~ 2 (log(n)+ < K >) -1 ), then the vector (7i,...,7 n ) is the needed 

approximation of the gradient. 

Notice that T{ = g§-Vk(^i, x n ) = ( n \y V(A, B, B), where the convex sets A = K \ 
and B = J2 i<i<n eVl Ky The randomized algorithm from [7J approximates V(A,B, ...,B) 
with the complexity 0(rj 4+0 ( 1 )e~( 2+0 ( 1 ) log(<5). This gives the needed approximation of 
the gradient with the complexity nO(n 8+ °^ (log(n)+ < K >) 2+0 W log((5 -1 )). 

3. Controlling the failure probability 5 . We need to approximate 0(n 2 log(Var(/))) 
values and gradients. To achieve a probability of success § we need that 



6 m j (n 2 (ni (log(n)+ < K >))) \ 



4 

This gives log((5) _1 ) ps 0(log(n) + log(log(n)+ < K >)). 



Remark 3.5: Let g(y\) = \og(VK.(e yi , ...,e Vn )) and g(y) = g(y) + h(y), \h(y)\ < a. We present 
here an alternative elementary way to approximate the partial derivative g'(yi). 
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Recall that the function g(.) is convex and < g"(x) < n : x G R. It follows that 

I g g{y)\<-o + — ( 25 ) 

The optimal value in (|25|) . <5 op i = 2..M, gives the bound 



i g(y + Sppt) - g{y) u\\s->i — fna\ 
I 7 g (y)\ < 2Vna. (26) 

Oopt 

The simple "estimator" (j25j) can be used instead of the interpolational algorithm from [7], but 
its worst-case complexity seems to be higher than that from [7]. | 



Theorem 3.6: Given a n-tuple K of well-presented convex compact sets in R n there is a 
poly-time algorithm which computes the number AV{K) such that 

Pr0b{1 ~ TfKT " 2 IJ Hh^Mf(i))) <2— } > .75 

l<i<ra 

T/ie complexity of the algorithm, neglecting the log terms, is bounded by 
0(n 12 (log(n)+ < K >) 2 ). 



Next, we focus on the case of Newton polytopes, in other words, polytopes with integer vertices. 
I.e. we will consider the mixed volumes V(P) = V(P±, P n ), where 

Pi = Hull ({vij : 1 < j < m(i), v itj G Z™}) . 



We define 

d(i) = min{A; : P,, C kHull(0, e\, e n )}, 

i.e. d(i) is the maximum coordinate sum attained on Pj. It follows from the monotonicity of 
the mixed volume that V(P\, P n ) < rii<j<fc^(*)- Such polytopes are well-presented if, for 
instance, they are given as a list of poly(n) vertices. This case corresponds to a system of sparse 
polynomial equations. Notice that the value V(P±, P n ) is either zero or an integer (BKK 
Theorem) and the capacity Cap(Vp) < ^ Y\i<i<k (inequality ([3])). 

The next theorem is proved in the same way as Theorem (|3.6p . 



Theorem 3.7: Given n-tuple o/P = (Pi, P n ) of well-presented integer polytopes in R" there 
is a poly-time algorithm which computes the number AV(P) such that 

Prob{l<—^<2 J] A(i,min(i,a//(P i )))<2— } > .75 
The complexity of the algorithm, neglecting the log terms, is bounded by 0(n 9 (n+log(Y\ 1<i<n di)) 2 ). 
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4 Open Problems 



1. Prove that for "random" convex sets V ^cap(v^)^ — nF^W w ^h high probability. This is 
true for the permanents of random matrices with nonnegative entries. 

2. The most important question is whether or not there exists a FPRAS algorithm for the 
mixed volume (or for the mixed discriminant). We conjecture that the answer is negative. 

3. Another important open problem is whether or not our mixed volume generalization (jlOp 
of Schrijver's lower bound on the number of perfect matchings in regular bipartite graphs 
[27] . [29] . [TT] is asymptotically sharp. 
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A Proofs of Theorems ( jZZD ) and (E3D 
A.l Useful (and well known) facts 

Fact A.l: Let it £ S n be a permutation and K = (K\, . . . , K n ) be a n-tuple of convex compact 
sets in R n . Then the next identity holds: 

V{K U ...,K n ) = V(K Al) ,...,K An) ) (27) 

Fact A. 2: We recall here the fundamental Alexandrov-Fenchel inequalities for the mixed vol- 
ume of n convex sets in R n : 

V(K u K 2 ,K^...,K n f > V(Ki,Ki,K 3 , K n )V(K 2 , K 2 , K 3 , ...,K n ) (28) 

Fact A. 3: Let {K%, .., Ki), i < n — 1; S, T be convex compact sets in R n . Define 

a Q = V{K 1 ,..,K i ,S,S,...,S),a 1 = V{K 1 ,..,K i ,S,S,...,S,T\a n ^ i = V{K 1 ,..,K i ,T,T,...,Ty 

Then the univariate polynomial U defined by U(t) = V(Ki, .., Ki, S + tT, S + tT, S + tT) is 
expressed as 

U(t)= (" ■ ')»!>■'■ (29) 

0<j<n~i \ J ) 

It follows from Facts (jA.ip and the Alexandrov-Fenchel inequalities that this univariate poly- 
nomial U is (n — i)- Newton. 

Fact A. 4: Let K = {K\, . . . , K n ) be a n-tuple of convex compact sets in R n . For a nonnegative 
vector (xi, ...,x n ), define the convex compact subset K = J2i<i<n x i^-i- Then the following 
identity holds: 

° -Vk(xi,...,x n ) = . 1 V(Ki, ...,Ki,K, ...,K) . (30) 



dx\...dxi ' ' {n — i)\ 
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Fact A. 5: If A is n x n real matrix and Aj > 0, 1 < i < n then 

V(\ l (AK l ),...,\ n (AK n )) = \det(A)\( H \i)V{K u . . . ,K n ) . (31) 

\<i<n 

Fact A. 6: . Let S, T; K2, K n be convex compact sets in R n . The next (additivity) identity 
holds: 

V(S + T, K 2 , K n ) = V(S, K 2 , K n ) + V(T, K 2 , K n ) (32) 

Fact A. 7: Let K = (K\ , . . . , K n ) be a n-tuple of convex compact sets in R n . Then the degree of 
the variable X{ in the polynomial Vk is degv K (i) = aff(i), where aff(i) is the affine dimension 
of Ki, 1 < i < n. 

A. 2 Auxiliary univariate inequality 
Lemma A. 8: 

Let R(t) = J2o<i<k a it l be an n — Newton polynomial, n > k and 

-1 



A(n,fc) = (min( 



SVnMX 



X 



where the polynomial sv ntk (x) = 1 + Ei<i<jfe(f )*(")• 
Then the following inequality holds: 



01 = R (0) > A(n, k) inf ^ . (33) 



Equality in {S^ is attained if and only if R(t) = R(0) (l + Ei<i<fc(f YQ) 
If n = k then 

\(n,n) = (mm = (^)" ' = : 9(n). (34) 

Equality in (3^> is attained iff R(t) = R(Q) (l + f ) n . 
Proof: Note that A(n, fc) < 1, k > 1. If R(0) = then clearly 

R '(0)=M^->X(n,k)M^-. 

Therefore we need to consider only the case R(0) > 0. Assume WLOG that R(0) = 1. It 
follows directly from the Newton inequalities ([8]) that 

a % < (" J ) .1 < / £ />. 
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Thus 

R(x)= aiX l <l+ (—) \ .) = sv(aix),x>0. 

0<i<k Ki<k ^ 71 ' W 



, c R{t) . fsv nk (a 1 x)\ . (sv nk {x) 

mi < mm = a\ mm ' 



It also follows that 



t>o t x>0 \ x J x>o V x 

which gives that 

R'(0) = ai > A(n,fc)mf— • 
The remaining statements can be now easily verified. I 

Corollary A. 9: Denote as Hom + (n,n) a convex cone of homogeneous polynomials with non- 
negative coefficients of degree n in n variables. Consider polynomials 

d 

p G Hom + (n,n); q G Hom + (n - 1, n - 1), q(x 1 , ...,x n -i) = - — p(x 1} ...,x n -i, 0). 

ox n 

Suppose that for all positive vectors X = (x%, x n _i) the univariate polynomials Rx{t) = 
p(x\, x n -i, t) are n-Newton. Then the following inequality holds: 

Cap(q) > X(n,deg p (n)) Cap(p) . (35) 

Proof: Note that because the coefficients of p are nonnegative the (univariate) degree 
deg(Rx) = deg p (n) for all positive vectors X G It follows from the definition of the 

polynomial q G Hom+(n — 1, n — 1) that 

q(xi, ...,x n -i) = R' x (0). 

It follows from the definition ([T]) of the capacity that 

Rx(t) = p(xi, ...,x n -i,t) > Cap(p)( Yl x i) t - 

l<»<n-l 

It now follows from the inequality (|33|) that 



x% , 

Ki<n-1 



q(xi, x n -i) = R' x (0) > A (n, deg p (n)) Cap(p) Yl 
which gives that Cap(q) > A (n, deg p (n)) Cap(p). I 
A.3 Proof of Theorem (E3D 

Proof: Let K = [K\, . . . , K n ) be a n-tuple of convex compact sets in R n . 
We associate with the Minkowski polynomial 

Vk G Hom + (n,n), V K (xi, ...,x n ) = V n {x\K\ + ... + x n K n ) 
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the following sequence of polynomials: 

dxi-\-\...dx 



q n = V K ;qi S Hom + (i,i),qi(xi, ...,Xi) = a f" \ q n (xi, x i} 0, 0), 1 < i < n - 1. 



Note that 

gi(x) = x dx ^ dx V K (0,...,0) = V(K 1 ,...,K n )x, 

and 

...,Xi) = ^ % + i(a?i,...,a;i,0),l <i<n-l. 

OXi+l 

Note also the obvious (but useful) inequality 

deg qi (i) < min(i, des 5n («)) = min(i, aff(i)) . (36) 

It follows from Fact (|A.4j) and Fact ()A.3j) that the polynomials qi,2 < i < n satisfy the conditions 
of Corollary (|A.9|) . Therefore 

Cap(qi) > A (i + l,deg qi+1 (i + 1)) Cap(q i+1 ) > A (i + l,min(i + l,aff(K i+1 )) Cap(q i+1 ). 

(37) 

(We use here the inequality (f36|) and the fact that A(n, fe) is strictly decreasing in both 
variables.) 

Multiplying the inequalities (|37|). we get the inequality ©. 



A.4 Proof of Theorem ( |2Jl 

Proof: (Proof of inequality ([3|) ), 

Sinse 



/i — 1\ 

A(i,fe) > A(i,t) = =: f — — j , 1 < k < i 
hence we get from Q that 

V(K 1 ,...,K n )>Cap(V K ) J] g(i) = ^Cap(V K ). 



2<i<n 



Proof: (Proof of the uniqueness part of Theorem (12. 

As remarked above, we follow in the present paper the proof of uniqueness in |17j . 
1. Assume that Cap(Vk) > 0. Suppose that I = aff(n) < n. As g(l) > g(n) hence 
V(K U ...,K n )> Cap(V K )g(l) J[ g(i) > ^Cap(V K ). 

2<i<n-l H 

Therefore if 

V(K 1 ,...,K n ) = —Cap(V K ) 
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then aff(n) = n. Using the permutation invariance (|27p . we get that 

a ff{i) = n,l < i < n. 

In other words the convex compact sets Ki all have nonempty interior. This fact together 
and the monotonicity of the mixed volume imply that all coefficients in the Minkowski 
polynomial Vk are strictly positive. 

2. Scaling. 

All coefficients in the Minkowski polynomial Vk are strictly positive, hence there exists 
an unique positive vector (a±, a n ) such that the scaled polynomial p = Vi aiKlj jllnK \ 
is doubly stochastic (see [TT]): 

^-KM,-,l) = l,l<*<n. 

We will deal, without loss of generality, only with this doubly stochastic case. 

3. Brunn-Minkowski. 

Let (z\, z n _i) be the unique minimizer of the problem 

min -,a?n-l)- 
x i >0,l<i<n-l;Jl 1 < i <„_ 1 Xi=l 

Such an unique minimizer exists as all the coefficients of q n -\ are positive. It follows from 
Lemma (|A.8|) and the proof of Lemma( |A.9|) . that 

VK(zi,...,z n - 1 ,t) = V(S + tK n ) = {at + b) n ,S= Yl z ^ 

\<i<n-\ 

for some positive numbers a, b. It follows from the equality case of the Brunn-Minkowski 
inequality [5] that 

K n = aS + {T n },a>0,T n eR n . 

In other words, K n = J2i<j<n-i A{n, j)Kj + {T n }, where A n j >0, 1 < j < n — 1 and 
V„ C /)'". 

In the same way, we get that there exist a n x n matrix A, with the zero diagonal and 
positive off-diagonal part, and vectors T\, ...,T n S R n such that 

K^Y^AihfiKj + iTi}. 

It follows from the doubly-stochasticity of the polynomial Vk that all row sums of the 
matrix A are equal to one. Indeed, using the identity (i30l) . we get that 

(n-1)! = (n-1)!— Vk(l, 1, 1) = V(SUM, SUM, SUM, Ki); SUM = Kx+...+K n , 1 < i < n. 

OXi 

As Ki = J2j^iA(hj) K j + i T i}i we g et ) usin g Fact(LA6J) and FactflAJiJ), that 
(n-1)! = V(SUM,SUM,...,SUM,^A(i,j)K j +{T i }) = 
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= ]T A(i,j)V(SUM, SUM, SUM, Kj) = (n - 1)! ]T A(i,j). 

Therefore J2j^iMhJ) = 1, 1 < « < n. 
4. Associate with the convex compact set K{ C R n its support function 

7i(X) = max < X,Y >,X G i? n . 

We get that 
As the kernel 

Ker{I - A) = {Y G i? n : (/ - A)Y = 0} = {c(l, 1, 1), c G i?}, 
it follows finally that 

ji(X) = a(X)+ < X,Lj >,X G R n 

for some functional a(X) and vectors L\,...,L n G i? n . This means, in the doubly- 
stochastic case, that Ki = K\ + {L^ — L{\, 2 < i < n. 



B Inequalities for Minkowski and Minkowski-like Polynomials 

Let / : — > be a differentiable positive- valued functional defined on the strictly positive 
orthant We assume that / is w-homogeneous, i.e. that f(ax\,...,ax n ) = a n f(xi,...,x n ) 

and the partial derivatives -£^f(x±, ...,x n ) > : (xi, ...,x n ) G -R++. We denote the set of such 
homogeneous functionals as PoH(n). 

We define the capacity as 

Cap(f) = inf J±h^l = inf f( Xl ,...,x n ). 

We define two subsets of PoH(n): 

Cav(n) - consisting of / G PoH{n) such that f~ is concave on all half-lines {X + tY : t > 
0} : X, Y G ; Vex(n) - consisting of / G PoH(n) such that /~ is convex on all half-lines 

{X + tY : t > 0} : X,Y G 

Recall theBrunn-Minkowski theorem : the Minkowski polynomial Vk(xi, ...,x n ) belongs to 
Cav(n). Therefore the results in this Appendix apply to the Minkowski polynomials. 

We also define the following Generalized Sinkhorn Scaling : 

SH(xi x)-(m y)-y- U?ll^2<± - X A x _i]k!^^A 

on \x\, x n j — \y\, yn) ■ Ui — a . — > ft — n, \ 

&:f(xi,...,x n ) li f{x 1 ,...,x n ) 
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Theorem B.l: If f £ Cav(n) then the following inequality holds: 

f(SH(x 1 ,...,x n )) < f(xi,...,x n ) . (38) 
If f £ Vex(n) then the reverse inequality holds: 

f(SH(xi, ...,x n )) > f(xi, ...,x n ) . (39) 

Proof: Let X = (xi,...,x n ) £ and Y = SH(xi, x n ). We can assume without loss of 
generality that f(X) = 1. If / £ Cav{n) then the univariate function g{t) = (f(X + tY))~ is 
concave for t > 0. Therefore 



9i t) < i9 (o) + = a + ^t r . 



n n 
We get, by elementary calculus, that 

d 

l<i<n 

The functional / is n-homogeneous, hence g{t) = t n f(Y + t~ 1 X) < (1 + t) n , and finally 
/(y + t- x X) < l^) n . Taking the limit t -> oo we get f(SH( Xl ,...,x n )) < 1 = /(s ls x„). 

The convex case is proven in the very same way. I 

Theorem (IB.lj) suggests the following algorithm to approximate Cap(f): 



X n+1 = Nor(SH(X n )) : Nor( Xl ,...,x n ) = (—,... ,^), a = J( TT n). (40) 

a a v l<*^« 

Corollary B.2: Consider f £ Cav(n). Suppose that Cap( f) > 0, 

log(Cap(/)) < log(/(x!, x n )) < Iog(Cap(/)) + e, < e < -L 

and Ili<i<n = 1; > 0, 1 < i < n. T/ien 

E (1 " "if'T,''' ) 2 S 10, • (41) 

l<i<n J{Xi,...,X n ) 

Proof: Let 7j = — f( Xl ,... x ) — • ^ follows from the Euler's identity that Xa<i<n7« = n an< ^ 
thus log(]li<j<n7i) < 0. 

Inequality (|38|) can be rewritten as 

/(— , — ) < f(x h ...,x n ). 

71 7n 
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Therefore \og{Cap(f)) < Cap(f) + e + log(rii<i<n 7i)i which gives the inequality 

-e < io g ( n t<) < o. 

l<i<n 

Finally, using Lemma 3.10 in [19] . we see that Xa<j<n(l ~~ 7*) 2 — 10 £ - I 

Corollary (|B.2p generalizes (with a much more transparent proof) corresponding results 
from H5 and [23]. 



T/ie following Lemma proves inequality Ittfy . 
Lemma B.3: 

1. Let p(t) = Eo<!<n a ^'i a i > be a polynomial with nonnegative coefficients. Assume that 
log(p(t)) is concave on -R++ and define q(x) = \og(p(e x )) . Then q(x) is convex on R and 
its second derivative satisfies the following inequality 

< q"(x) < n . (42) 

2. Let p(t) = J2o<i<n a it 1 -, a i > be a polynomial with nonnegative coefficients. Assume that 
p(t)m,m > n is concave on R + . Then 

< q"{x) < f(n,m), (43) 

2 

where f(n, m) = n — ®- ifn < y ; and f(n, m) = y otherwise. If n = m then f(n, m) = j 
and the upper bound |j5p is attained on polynomials p(t) = (a + tb) n ; a, b > 0. 

Proof: 

1. The convexity of q(x) is well known. The concavity of log(p(i)) is equivalent to the 
inequality (p'(t)) 2 > p(t)p"(t) : i > 0. Putting y = e x , we get that 



q "( x \ = p"(y}y 2 + p'(y)y _ ( p'(y)y \ ' 

p{y) p(y) \ p(y) J 



p(y) p{y) V p(y) 

The concavity of log(p(t)) gives that p ^ ( VfrffiO — ^' 

The coefficients of the polynomial p are nonnegative, hence P pffi < n - This last observa- 
tion proves that 

,, ( x . p'(y)y . 

q (x) < — — < n. 

p(y) 

2. Our proof of (|43p is a direct adaptation of the above proof of (|42p . We use the following 
characterization of the concavity of p(t)™,m > n: 
(P'(t)) 2 > ^iP{t)p"{t) : t > 0. 

I 

We note that just the nonnegativity of the coefficients implies the quadratic bound 

q"{x) < 0.25n 2 . 
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